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Abstract 

Quantum fluctuations of a scalar field and its derivatives are calculated when the 
field is confined between two parallel plates satisfying Dirichlet or Neumann boundary 
conditions. After regulation these fluctuations diverge in general when one approaches 
one of the plates. The energy density and the pressure between the plates is only 
consistent with the total Casimir energy when the canonical energy-momentum tensor 
is augmented by the Huggins term so to satisfy the requirement of conformal invariance 
for a massless, scalar field. 

m : 

in ; 1 Introduction 

Fluctuations in a quantum field can under certain conditions give rise to forces between 
macroscopic objects. This was first realized by Casimir who considered an electromag- 

^^O ■ netic field between two parallel plates separated by a distance L. From the quantum 

fluctuations of the field described by the Lagrangian Cem = (1/2)(E 2 — B 2 ), he found an 

Q ! attractive force corresponding to the potential energy 

O ' vr 2 

per unit plate area. Assuming a uniform distribution of the energy between the plates, 
£h ' one then has in this region the energy density 

3 : TT 2 

FT S ° = -72MA (2) 



while outside of the plates one expects it to be zero. 



That it should be so, is not entirely obvious. For the similar Casimir energy inside a 
spherical shell, explicit calculations show that the energy density varies with the distance 
from the shell |2j, both inside and outside. With two parallell plates, the energy density is 
given by the expectation value of the fluctuations in the electric and magnetic fields, 



(E 2 ) + (B 2 ) (3) 



These fluctuations were first calculated by Liitken and Ravndal 3 and found to vary 
between the plates. In fact, both of them diverge when one of the plates is approached. 
But this position dependence cancels in the sum, resulting in a constant energy density. 
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As pointed out by Brown and Maclay[I], this follows from the more fundamental require- 
ment that the energy-momentum tensor T^ u of the electromagnetic field is traceless since 
the photon is massless. With two parallel plates normal to the unit vector n M = (0; 0, 0, 1) 
along the z-axis, its expectation value must then have the form 

vr 2 

{T lxv )=- -^jj (flu* + 4? V^ ) ( 4 ) 

with the metric rjau = diag(l, — 1, — 1, — 1). By construction, the energy density is then 
given by (Too ) and the pressure on the plates is (T zz ) = — 7r 2 /240L 4 . 

A detailed calculation of these effects for the electromagnetic field is a bit cumbersome since 
it must be expanded in transverse electric and magnetic multipoles. This complication is 
absent for the massless scalar field (f>(x) which is spinless and described by the Lagrangian 
Cq = ^(d^(f>) 2 = \[4> 2 — (V0) 2 ]. In fact, it is usually the scalar field which is used when one 
wants to arrive at the Casimir energy in the simplest way. Not unexpected, one obtains 
a total energy which is exactly one half of the electromagnetic result © since the photon 
has two spin degrees of freedom. The corresponding energy density should then follow 
from the Hamiltonian density, 



p(o) _ I 
L ~ 2 



(0 2 ) + (W) 2 ) (5) 



These expectation values giving the field fluctuations between the plates will be calculated 
in the following. As for the electromagnetic system, both of the contributions to the energy 
density vary with the position between the plates. However, in this scalar case their sum 
does not add up to a constant value. Instead we find that the sum diverges near the 
plates. As a consequence, the integrated energy density does not agree with the known, 
full energy and is in fact divergent. 

The reason for this complication has been known for a long time and is caused by the 
conformal invariance which is present for all free, massless fields and requires a traceless 
energy-momentum tensor as for the electromagnetic field. In the above calculation we used 
the Hamiltonian density which is the 00-component of the canonical energy-momentum 
tensor 

T^l = d^d v <P-^(d x <P) 2 (6) 

which has the trace T^* 1 = — {d\4>) 2 . The improved energy-momentum tensor derived by 
Callan, Coleman and Jackiw[S] has an additional piece 

AT^ = -^8 U - V ^d 2 )^ 2 (7) 

which had been obtained earlier by Huggins 6 . The full tensor T^ v = T$ + AT^ U has 
zero trace using the equation of motion d 2 (j) = 0. We then find a resulting energy density 
which is constant and integrates up to the full energy. The need for this extra term was 
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first pointed out by deWitt who calculated the full energy-momentum tensor and showed 
that it gave a finite Casimir energy [7]. More recently, it has been calculated by Milton 8 
using Green's functions methods, obtaining the same result for the energy density. 

In the next chapter the scalar field between two parrallel plates is quantized in a stan- 
dard mode expansion for both Dirichlet and Neumann boundary conditions. The Casimir 
energy is calculated making use of combined zeta-function and dimensional regulariza- 
tions to obtain a finite, physical result. With these quantized modes one can then also 
calculate the fluctuations in the field and its derivatives. This is done in the following 
chapter where the same regularization suffice here to make such local quantities finite. 
We can then also calculate the energy density and it is shown that the conformally co- 
variant energy-momentum tensor gives a finite and consistent result both for Dirichlet 
and Neumann boundary conditions. In an appendix we give a short summary of confor- 
mal transformations and the derivation of the correct energy-momentum tensor from the 
invariant theory in curved spacetime. 



2 Eigenmodes and the Casimir energy 



The massless scalar field <j) = </>(x, t) is defined in a volume closed by two parallel infinite 
plates. It is convenient to introduce the coordinate split x = (x-p, z) where the x^ = (x, y) 
describes the position along the surfaces, while the coordinate z describes the position 
normal to the plates and therefore z e [0, L] in the region between them. The dynamics of 
the classical field is given by the Klein-Gordon field equation d 2 (p = 0. At the surfaces 
we will impose Dirichlet <j>\ z =o,L = or Neumann d z (f)\ z =o t L = boundary condition. We 
then have the following two sets of orthonormalized eigenmode functions 

«-<*•«> - £(££2j)*~-* (8> 

where 

w = (4 + *£) 1/2 

is the frequency of the mode and k n = nir/L with n = 1,2,3,... is the longitudinal 
wavenumber. The Neumann mode with n = is constant in space and doesn't contribute 
here. 



The field can now be expanded in terms of these basic modes 



oo „ 

(M) = £ / 

n=l 



d 2 k T fl 



(2vr) 2 V 2lo 



flk T .n'Uk T ,n( x > t) + a j CTjn - Uk T ,n( x ; 



0) 



where the expansion coefficients become creation and annihilation operators after quan- 
tization. With the above normalization of the eigenmodes, they satisfy the canonical 
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commutator 



[ak T ,n,a\^ n/ ] = (2vr) 5(k T - k^) 5 n ^ (10) 

We will here in particular consider the vacuum state 1 0) defined by the standard condition 
a k T ,n|0) = for all modes specified by the mode quantum numbers (kr,ra). 

Summing up the zero-point energies (l/2)u> of each mode, the full vacuum energy per unit 
plate area is given by 

^^£/Ss(^ + (Wi) 2 ) 1/2 (id 



n=l' 



In order to make this divergent expression finite, we make use of the basic integral 

ddk 1 _ nN-d/2) 2)d/2 _ N 



(2ir) d (k 2 + m 2 ) N {^) d / 2 T{N) 

from dimensional regularization. This is convergent when the power N is suffiently large 
compared with the dimension d. When this is not so as in the above integral for the 
vacuum energy, we define its value by analytical continuation in the dimension d. In our 
case d = 2 and N = —1/2 which gives 

E o = -J^sE n3 ( 13 ) 

n=l 

The remaining, divergent sum is done by analytical continuation of the Riemann zeta- 
function, defined by 

oo 



coo = £ t. 



1 n° 

n=l 



for Res > 1. This gives Y^=i n ^ = C( — 3) = 1/120 and the scalar Casimir energy takes 
on the well-known value 



T 2 



E ° = -I440P (14) 

corresponding to the pressure p = — 7r 2 /480-L 4 normal to the plates. We thus have the 
same result both for Dirichlet and Neumann boundary conditions. In other regularization 
schemes there sometimes arise a subtle difference between the two having to do with the 
zero mode n = which doesn't contribute in dimensional regularization [§]. 



3 Field fluctuations 



Since the field is linear in creation and annihilation operators, the expectation value of the 
field is zero in all states. This applies in particular to the vacuum state, (0\4>(x)\0) = 0. 
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But the fluctuation (0\4> 2 (x)\0) is non-zero. It can be calculated from the field expansion 
© which gives 

/ , ,i, m v 2 ^2, f f d 2 k T d 2 k' r 1 t 

<O|0V)|O) = I J2jj—^ — {0 \a kT „al T J ) 

x / sin(fc n ^)sin(A; n /2;) \ i(kT _ k y. XT _ i(aJ _^ )t 
\ cos(A; n z) cos(A; n /z) / 

1 ™ f d 2 k T 1 . 

n=l v ' 

where = ttz/L gives the distance from one plate. Here and in the following upper signs 
are for Dirichlet and lower signs are for Neumann boundary conditions. The momentum 
integral can now again be made finite with dimensional regularization as in the previous 
chapter. One then finds 

oo 

(0|^ 2 (x)|0) = -_£„(l =F cos 2n0) (15) 

n=l 

The first sum here is given by Y^n=i n = —1/12 with zeta- function regularization while 
the second sum becomes 



oo \ d 1 / d \ 1 

En cos 2n6 = - — sin 2n6 = — I — ) cot 9 = 
2 86^ A\d6 Asm 2 6 
n=i n=i y ' 

We thus have for the field fluctuation 

WWIO) = ^(1^) (16, 
It is seen to diverge near the plates where 6 — > 0, it. 

The field fluctuation outside just one plate can now be obatined by removing the other 
plate to infinity, i.e. letting L — > oo in the above result. We then find 

(^ 2 WI«)=T^ (17) 

It is difficult to imagine any physical measurement of these scalar fluctuations. But the 
corresponding electromagnetic field fluctuations will perturb the energy levels of any atom 
near such a plate and should in principle be observable 3 . 



4 Energy-momentum tensors 

From the canonical energy-momentum tensor Q we have the expression © for the energy 
density between the plates. The first expectation value ( <j) 2 ) can be calculated as in the 
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previous chapter. We then obtain 



-i oo „ j2j. 

<0|^|0) = 2lE J T^"(lTco S 2n9) 



n=l 



After dimensional integration over the transverse momenta, which is the same as for the 
full Casimir energy we have 



TT 2 °° 



(O|^|O) = -^ ZI ^n 3 (l T cos2n0) 



71=1 



The remaining sum is similarly found to be 



00 1 / 8 \ 3 1 

£ n 3 cos2n0 = --(-) cot*=-/(0) 

n=l ^ ' 



where the function 



/(*) = — a ~— ft (18) 
sin sin ft 

is the same as appears in the derivation of the electromagnetic fluctuations in the same 
geometry jS]. We thus have (0|^ 2 | 0) = —(A^pB) with 

^ = T^i £- 3 = 7440ZI ^ 

n=l 



and 

TT 2 °° TT 2 

B = wji E n3 cos = 9^4 / W (2°) 

71=1 



The other needed expectation values can be obtained in the same calculational scheme with 
the result (0|(3^) 2 |0) = -3(A±B) and (O|(V T 0) 2 |O) = 2{A^B). Thus (0 1 (<9 A (/>) 2 | 0) = 
±6B. 

From (JSJ) we now find £0 = — (^4 ± 21?) for the canonical energy density. Due to the 
presence of the 5-term, it is not constant between the plates and actually diverges when 
one approaches one of the plates. In the same way we will find that the expectation value 
( T zz ) does not reproduce the force between the plates. However, let us now consider the 
Huggins term (|7j). Using the equation of motion d 2 4> = 0, it can be written as 

AT^ = —(d^tff + <HA</> - V^(d x <P) 2 ) (21) 

Since ( 4> 2 ) = —((ftd 2 (j)), the first two parts cancel for the 00-component which thus be- 
comes ( AToo ) = ±2-B and cancels the position-dependent part of £q. The full, conformally 
covariant energy density is therefore simply £ = —A as expected from the full energy (|14|) . 
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The pressure follows similarly from zz-component of the full energy-momentum tensor 
which can be simplified to 



Notice that ( (<9 z </>) 2 ) cannot be directly related to (4>d^(p} = 3(A =p B). Putting in the 
other expectation values, we find (T zz ) = —3A = — 7r 2 /480L 4 which is the presure on the 
plates. 

Needless to say, the conformal Huggins term will also contribute to the energy density 
and pressure when the system is in thermal equilibrium at finite temperature. The above 
vacuum expectation values are then replaced by thermal averages. These can also be 
calculated by standard methods to give results on closed form since the fields are free, but 
this will not be pursued here. 

5 Conclusion 

For the scalar contribution to the Casimir force of two parallel plates the Huggins term is 
not important since it integrates out to zero in the total energy. But for the local energy 
it is essential for consistency. We have seen that the quantum fluctuations vary between 
the plates and actually diverge when one of the plates are approached. It is only with the 
conformal Huggins term that the energy density is evenly distributed between the plates 
and zero outside. 

Scalar fields play an important role in modern cosmology|llj. At very early times the 
inflationary epoch is driven by the vacuum expectation value of the scalar inflaton field. 
And at much later time one can model the dark energy responsible for the acceleration 
of the cosmic expansion by a scalar quintessence field. In both of these cases there is 
the question about what the kinetic energy of the field should be. Quantum effects for a 
scalar field in a curved background will generate a counterterm term Rcj) 2 which should 
therefore be included in the Lagrangian. It will generate a term in the energy-momentum 
tensor proportional to the Huggins term although there is a priori no underlying conformal 
invariance in the system. For this reason, such a term should in general be included but 
with an unknown coefficient. It will modify the behaviour during the inflationary epoch [T2*] 
and later affect the acceleration of the universe. In fact, recently is has been shown that 
an effective equation of state p = wp for the dark energy can result with w < — 1 which 
represents very strange physics now called 'phantom energy' [T3] . A better understanding 
of these effects will have to wait until we have a deeper theory of fundamental interactions 
where such scalar fields arise naturally. 




(22) 
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6 Appendix 



Conformal invariance and the energy momentum tensor for scalar fields is discussed in the 
book by Birrell and Davis 10]. We will here give a short summary and a bit more detailed 
derivation of the Huggins term. 

For a system described by the Lagrangian £, the action is 




in curved spacetime with the metric g^ v with g = detg^ u . Under an arbitrary variation 
9^u -> 9fj,u + Sg^v it changes by 

5S = j 'd 4 x((J/=g£ + </=s<J£) 
Since 5^/—g = —\^J —gg i iv&g^ v we can write the change in the action as 

iT^SgT (23) 

- g„uC (24) 
is the energy-momentum tensor. For the minimal Lagrangian 

£ (0) = ynvQ^dvcf) this 
gives the canonical energy-momentum tensor © in flat spacetime. 

Under a conformal transformation g^ — ► Q 2 g^ where Q = Q(x), the metric changes by 
an amount 5gau g^v If we now demand invariance under such a transformation, we see 
from (j2IIJ) that 5S = then requires = 0, i.e. the energy-momentum tensor must be 
traceless. 

Since the scalar field has mass-dimension cL = 1 in a 4-dimensional spacetime, it will 
transform as <j> — ► $7 _1 <^. In the kinetic energy we will therefore have the change 

where the comma derivative <j> „ = du<i>- This extra piece generates new terms in the 
kinetic energy of the field which must be cancelled by an additional part in the minimal 
Lagrangian. We will now show that the extended Lagrangian 

C^g^d^d^+^m 2 (25) 

where R = g^ v R^ u is the scalar Ricci curvature, is conformally invariant. Since it follows 
directly from the definition that yj— g — ► Q\j—g under the transformation while for the 
curvature 10 

r -> n~ 2 R - 6fi~ Vfy/«/ , 



where 



SS = | / d 4 x, 



T — 2 

" dg^ 
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we see that the full Lagrangian changes into 
The extra term is 



-gQ 



Q 



1 



[IV 



after having changed all comma derivatives of scalar quantities into covariant derivatives. 
Now we can simplify this expression to the total derivative Q = — \{Q.~ 1 Q,'^(j) 2 )^ which 
gives a vanishing contribution in the integral for the action. The Lagrangian (|25|) is 
therefore invariant under conformal transformations. 

Under the variation — * <7„„ + 8g^ v we see that the last term in the invariant Lagrangian 
(|25[) will generate a new contribution to the energy-momentum tensor from the additional 
term 

5SR = ~kl ^ x R + ^9 SgTR^ + V=5 gTSR^ 

in the variation of the action. Here we will need the variation of the Ricci tensor which 
can be expressed in terms of the Christoffel symbols as 



R 



[IV 



/iv,a 



■pet _j_ "pQf t~*P "p/3 "po; 

fia,v fii> af3 fJ.cn v(3 



(26) 



The calculation is most easily performed in Gaussian normal coordinates, where the 
Christoffel symbols are zero, but their derivatives are non-zero. We then have 



gTSR 



[IV 



(27) 



[iv,a [ia,vJ \i> ui [iv i> u [ivl ,ct 

since the derivatives of the metric is zero in these coordinates. The expression within the 
parenthesis is a vector w a since = g a/3 (6g j 3 fl)U — \8g^ v ^) is a tensor. Together with 
ST" = \g v ^5g v p^ we then have for this vector 

w a = -g^g^[Sg^-6g p ^} 
Collecting terms, we can write the new contribution to the variation as 

8SR = ~kS dV^^V" - g^g^^g^ - Sg^, av . 

where E^ v = R^ v — ^g^R is the Einstein tensor. In the last term we perform two partial 
integrations resulting in 



SS R = ^J d A x^g-(E^<p 2 - [(</>V - g^ 2 y a a ]yg 



[IV 



Comparing with the definition 1)23(1 of the energy-momentum tensor, we see that we now 
can identify the general Huggins term. In a general coordinate system where partial 
derivatives are replaced with covariant derivatives V^, it is 



AT, 



[IV 



1 

6" 



-E, 



[iv\ 



y^v 



g^v 



□ 



(28) 



9 



where □ = VVj, is the d'Alembertian operator. In Minkowski spacetime it simplifies to 
Q) in the text. 

The full energy-momentum tensor 

is now traceless, 

n = {d^f -2(d^f + 1 -{R-2R)4> 2 - -4n)0 2 

= -{d^f -^R<P 2 + {d^f + = 

since □ = (l/6)R(fi is the equation of motion for the conformally invariant Lagrangian 
(|25jl. An additional massterm oc 4> 2 breaks the invariance, while a oc 4 coupling preserves 
it and thus will give rise to a traceless energy-momentum tensor for the interacting field. 



References 

[1] H.B.G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948). 

[2] K. Olaussen and F. Ravndal, Nucl. Phys. B192 237, (1981). 

[3] C.A. Liitken and F. Ravndal, Phys. Rev. A31, 2082 (1985). 

[4] L.S. Brown and G.J. Maclay, Phys. Rev. 184, 1272 (1969). 

[5] C.G. Callan, S. Coleman and R. Jackiw, Ann. Phys. N.Y. 59, 42 (1972). 

[6] E. Huggins, Ph.D. thesis, Caltech, 1962 (unpublished). 

[7] B. deWitt, Phys. Rep. 19C, 295 (1975). 

[8] K.A. Milton, Phys. Rev. D68, 065020 (2003). 

[9] F. Rsivndal, \hep-ph/000920Si 



[10] N.D. Birrell and P.C.W. Davies, Quantum fields in curved space, Cambridge Univer- 
sity Press, 1982. 

[11] T. Padmanabhan, Phys. Rep. 380, 235 (2003). 

[12] F. Lucchin, S. Matarrese and M.D. Pollock, Phys. Lett. B167, 163 (1986). 
[13] V. Faraoni, Phys. Rev. D68, 063508 (2003). 



10 



